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Abstract. Massive quarks are included in the Curci- Ferrari model and the theory is renormal- 
ized at two loops in the MS scheme in an arbitrary covariant gauge. 
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Recently there has been renewed interest in examining covariantly gauge fixed QCD where 
the gluon and ghost fields are given explicit mass terms. For instance, in |l], ^, || such models 
have been used to investigate how a mass gap emerges for the gluon. In particular a dimension 
two operator obtains a non-zero vacuum expectation value which generates a gluon mass. Indeed 
in this context the effective potential calculation of |J demonstrated that the non-perturbative 
vacuum favoured a non-zero vacuum expectation value for the simplest dimension two operator 
possible in Yang-Mills theories. In the main these studies did not involve quarks and were 
effectively based on a version of Yang-Mills theory with massive gluons originally introduced 
by Curci and Ferrari in Q. There a renormalizable theory of massive gluons was constructed 
with the aim of being an alternative to the Higgs mechanism for endowing vector bosons with 
mass. The main shortcoming of the Curci-Ferrari model, however, was the breaking of unitarity 
directly as a result of the massive gluon, ||, ||, 0, |||, ||, 1C]. Nevertheless the model has proved 
useful for a variety of reasons. From a field theoretic point of view it is of interest due to the 
nonlinear nature of the gauge fixing term which introduces a quartic ghost self-interaction as 
well as modifying the usual ghost gluon interaction present in a linear covariant gauge fixing, 

The nonlinear property has been examined in [10, 11|. Moreover, there has been a debate 
on the BRST symmetry of the Curci-Ferrari model. For instance, Ojima, ||, has carried out a 
comprehensive examination of the BRST algebra with and without a mass term and explicitly 
constructed a negative norm state which therefore supports the lack of unitarity in the model. 
Other such states have been determined in j|]. Another feature which emerged in these papers 
was the non-nilpotency of the BRST charge which appears to follow as a consequence of the exact 
form of the Lagrangian. A more recent study of this has been given in [12]. A final motivation 
for considering massive gluons rests in phenomenological considerations. For example, lattice 
results, [pi ) suggest that at low energies the gluon is massive though the precise form of the 
propagator at this scale is not known explicitly but is clearly dependent on non-perturbative 
properties. However, one can use a Curci-Ferrari type theory with the hope that it can provide 
useful insight into physics which is in someway dependent on a gluon which is massive. For 
instance, diffractive scattering has been examined in [14] with such a motivation. 

Whilst the majority of the papers in this area concentrate on Yang-Mills theory with a 
massive gluon the real world is based on QCD which involves quarks. Therefore, if complete 
studies of the property of a massive gluon are to be performed, massive quarks need to be 
included. In Q matter fields were considered but were shown not to affect the failings of 
the Curci-Ferrari model such as lack of unitarity or non-nilpotency of the BRST charge. The 
multiplicative renormalizabilty of the model, which was proved in |||, ||, || [I(J, was unaffected. 



Given that the Curci-Ferrari model has been renormalized at one loop in M, 12, uM and more 



recently at two loops in MS in [16], the purpose of this letter is to extend the renormalization 
of the Curci-Ferrari model at two loops to the case where massive quarks are included. This is 
important for various reasons. First, the main motivation of [16] was to provide a calculational 
tool for renormalizing Green's functions of Yang-Mills theories where one could not naively 
nullify the momenta of several external legs without introducing spurious infrared infinities. 
Ordinarily one can handle such infrared problems if some form of infrared rearrangement, |L"5| , 
16], is performed to eliminate those divergences which cannot be distinguished from ultraviolet 



ones in a dimensionally regularized calculation. Moreover, given that infrared rearrangement is 
usually carried out by hand it does not lend itself readily to automatic calculations by computer. 
One simple resolution of the infrared problem is to introduce an infrared mass regularization. 
However, in using the Curci-Ferrari model, Q, which naturally incorporates a gluon and ghost 
mass and preserves renormalizability one resolves the infrared problem, [|l9| , |l0| , [la], and opens 
the path to automatic calculations. Therefore, extending the Curci-Ferrari model to include 
quarks is the natural way to proceed in order to provide a tool similar to [14] for QCD. However, 
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there is an additional practical point of view for our work. In conventional QCD perturbative 
calculations the fields of the Lagrangian are massless. However, in the real world the quarks have 
a physical mass. Therefore, including an explicit mass term for the quarks which is independent 
of the gluon mass we will have a Lagrangian which in principle includes as a special case the 
more realistic situation of massive quarks and massless gluons and ghosts as well as allowing 
us to interpolate between various different scenarios. For the two loop MS renormalization 
of the Lagrangian we do not expect significant differences with conventional results for the 
renormalization group functions. However, the tool we provide here will, for example, be useful 
in studying the mixing of operators under renormalization to operators of the same and lower 
dimension since we will have an explicit quark mass at our disposal. 

We take as our Lagrangian, H, f|. 
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where 1 < a < Na, 1 < I < Np with Np and Na the dimensions of the colour group fundamental 
and adjoint representations respectively where the structure constants are f abc , 1 < i < Nf where 
Nf is the number of quark flavours, g is the coupling constant, m is the gluon mass and hence 



the basic mass scale of the classical theory and c a c b 
G^, follows from the definition of the covariant derivative as Gr" 



c a df,c b 



d^Al 
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We have included the usual covariant gauge fixing term with parameter a which, to ensure that 
the action is BRST invariant, [Q, [6|, ||, dictates the form of the ghost interactions with 
a gluon of mass m and a ghost of mass y/am. To make the two loop calculations easier to 
perform we have chosen to parameterize the quark mass with the parameter j3 so that there is 
one basic mass parameter. This means that (3 will get renormalized but the full quark mass, 
y//3m, will be renormalized in such a way that it is independent of a. The case of massive 
quarks but massless gluons is recovered by setting m — > in such a way that \f]3m remains 
finite. With this Lagrangian the Minkowski space propagators for the gluon, ghost and quark 
fields are respectively, 
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The renormalization of (||) proceeds along usual grounds. First, we introduce the renormal- 
ized variables via 

A^ = JZ~ A A a » , c% = ^JY c c a , c£ = v^c" , = 
g Q = Z g g , to = Z m m , a = Z~ x Z A a , (3 Q = ZpP (3) 

where the subscript, > denotes bare quantities. Since our calculation builds on the Yang- Mills 
version of the Curci-Ferrari model we have not assumed the usual Slavnov- Taylor identity of 
Z a = 1. Therefore, we have seven independent renormalization constants to compute. However, 
we do have various cross-checks on the results we will obtain for the renormalization constants. 
Given that the parameter m appears in both the gluon and ghost sectors we ought to obtain 
a consistent renormalization for it from considering independently the gluon and ghost two 
point functions. Moreover, we can check that our wave function renormalizations are correct by 
examining the various vertex corrections. The same coupling constant renormalization constant 
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ought to emerge in all cases. Given the explicit values of the renormalization constants the 
corresponding renormalization group functions follow from the relations 
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where 7/3(0) is defined by 
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and a = g 2 /(I6ir 2 ). (The parameter (3 ought not to be confused with the /3-function, /3(a).) In 
ordinary QCD one would have 7a(o) = — 7a (a) but since Z a ^ 1 in the Curci- Ferrari model we 
do not expect this to be restored in the presence of quarks. Moreover, this is the origin of the 
second terms in the expressions for 7 m (a) and 7/3(0) which will be dependent on a. 

Given that we are renormalizing a massive version of QCD we need to have an algorithm 
for computing the ultraviolet divergences of the massive multiscale two loop Feynman integrals 
which contribute. First, we use dimensional regularization with d = 4 — 2e. Second, we follow 



the strategy of [ 2C , 21 ] where massive Feynman integrals are expanded in powers of the external 
momenta based on the identity 
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The expansion is terminated by the rule that when the powers of momenta exceed those which 
can appear in the Green's function through renormalizability, then they are dropped. So, for 
example, the gluon propagator is only expanded to 0(p 2 ) with the 0(p 3 ) terms being dropped 
where p is the external momentum. Consequently, one is left with massive vacuum bubble graphs 
where because of the different masses of (|]), not all the masses are equal. However, massive 
vacuum two loop bubbles have been evaluated to the finite part in p2^| , for example, though 
we only require the Laurent expansion to the simple pole in e. For instance, the basic vacuum 
bubble with three different scales is given by 
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where 7 is the Euler-Mascheroni constant and J k = J d d k/(2it) d . The expression for different 
powers of the propagators are determined by differentiating with respect to the parameter a, (3 
and m 2 . It is worth noting that the appearance of In a and In (3 terms could in principle lead to 
a non-analytic renormalization constant. However, these ought to cancel when all contributions 
from one and two loop graphs are included. In addition, in reducing the integrals to vacuum 
bubbles, partial fractions have been used which can give rise to other potentially singular terms. 
For instance, 
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This provides another internal check since there are no singular terms in the original Lagrangian 
and for (||) to be renormalizable they ought not to remain after the two loop calculation. Given 
the large amount of algebra which arises due to the expansion to vacuum bubbles, to handle 
their evaluation for different masses we have written an algorithm in a symbolic manipulation 
language, Form version 3, p^ ]. The calculation proceeds automatically, since, for example the 
Feynman diagrams are generated using the Qgraf package, [24]. Moreover, the renormalization 
constants are extracted by computing the Green's functions in terms of bare parameters then 
rescaling by (||) after the pole structure has been determined and following the procedure of 
[25]. This reproduces the usual method of subtractions automatically. 

We have renormalized the gluon, ghost and quark two-point functions as well as the three 
3-point vertices. We obtain the following renormalization constants 
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where T a T a = C F , f acd f bcd = C A 5 ab and Tr (T a T b ) = T F 5 ab . We have recorded these explicitly 
since the double pole in e follows from the form of the one loop simple pole and therefore 
provides another check on our computation. Moreover, they are /3-independent since we use 
a mass independent renormalization scheme. Prom these values we obtain the renormalization 
group functions at two loops in MS, 
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The expression for the (3- function agrees with the scheme independent results of [26, [27j]. The 



renormalization group functions for the wave functions of the fields and a agree with the corre- 
sponding results of [ 



26, 27 



ja(o) + 7 Q (a) = a 



in the Landau gauge, a = 0. Indeed 
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which is independent of Nf and vanishes when a = 0. Moreover, ([ll]) is equivalent to the 



statement that the ghost gluon vertex does not get renormalized in the Landau gauge, [29 



which was originally verified at one loop for © in [§9|. With the presence of quarks the 
gluon mass dimension now depends on Nf at one loop. However, to check we have correctly 
determined the quark mass anomalous dimension we need to compute the anomalous dimension 
of m„ = \fj5m. From 
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which agrees with the two loop MS result of [31, 32] in our conventions and is independent of 



a. For completeness, the ghost mass dimension is determined from 

m 2 = am 2 (14) 

which implies 

7m c 0) = 7m(«) + \la{a) (15) 

giving 
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which becomes Nf and gauge dependent at two loops. Further checks on the correctness of 
7 m (a) are that the Yang- Mills sector of our expression agrees with the three loop result of ||] 
and the one loop expression of |[3| where the multiplicative renormalizability of the composite 
operator ^A^A® 11 — ac a c a was verified by determining the mixing matrix of the renormalization 
of the constituent dimension two operators. For completeness, we have evaluated 7(a) in the 
Landau gauge for QCD, N c = 3 and Tp = i, and found 

2 

7 m (a)| Q=0 = (8^-105)^ + (548^-4041)^ + 0(a 3 ) . (17) 

It is worth noting that both coefficients are negative for Nf < 8 which implies that the gluon 
mass runs to zero in the ultraviolet limit in this case. 

In conclusion we have provided the full two loop renormalization of the Curci-Ferrari model 
with massive quarks in the MS scheme. Whilst the nature of the model may appear unphysical 
with a gluon mass and lack of unitarity, it is important to recall that one long term aim is 
to use the the model to attack the renormalization of other Green's functions where quark 
mass effects will become important and a natural infrared mass regularization is necessary. 
Further, it is worth noting that models with non-abelian symmetries and massive gluons may 



have important phenomenological consequences. For instance, it has been shown in [33] that 
including 0(l/Q 2 ) power corrections in the operator product expansion improves the matching 
of calculated physical quantities with experiment. Although such corrections derive from a gluon 
mass, albeit tachyonic in origin, they indicate the potential importance of massive gluons, such 
as that of (||), to probe non-perturbative physics. 

References. 

[1] K.-I. Kondo, Phys. Lett. B514 (2001), 335. 



[2] K.-I. Kondo, |hep-ph/0110015 . 



[3] H. Verschelde, K. Knecht, K. van Acoleyen & M. Vanderkelen, Phys. Lett. B516 (2001), 
307. 

[4] G. Curci & R. Ferrari, Nuovo Cim. A32 (1976), 151. 

[5] G. Curci & R. Ferrari, Nuovo Cim. A35 (1976), 1; G. Curci & R. Ferrari, Nuovo Cim. A35 
(1976), 273; Nuovo Cim. A47 (1978), 555. 

[6] I. Ojima, Z. Phys. C13 (1982), 173. 

[7] L. Baulieu, Phys. Rep. 13 (1985), 1; R. Delbourgo, S. Twisk & G. Thompson, Int. J. Mod. 
Phys. A3 (1988), 435; T. Hurth, Helv. Phys. Acta 70 (1997), 406. 



7 



[8] J. de Boer, K. Skenderis, P. van Nieuwenhuizen k A. Waldron, Phys. Lett. B367 (1996), 
175. 



[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17; 

[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 

[27; 

[28 

[29 
[30 
[31 
[32 

[33 
[34 



F. Delduc k S.P. Sorella, Phys. Lett. B231 (1989), 408. 
A. Blasi k N. Maggiore, Mod. Phys. Lett. All (1996), 1665. 
R. Delbourgo k PD. Jarvis, J. Phys. A15 (1982), 611. 
E. Andriyash k K. Stepanyantz, |hep-th/0204211| . 

Ph. Boucaud, A. Le Yaounac, J. P. Leroy, J. Micheli, O. Pene k J. Rodriguez-Quintero, 
Phys. Rev. D63 (2001), 114003. 

J.R. Forshaw, J. Papavassiliou k C. Parrinello, Phys. Rev. D59 (1999), 074008. 



V. Periwal, |hep-th/95090&l . 

J. A. Gracey, Phys. Lett. B525 (2002), 89. 

A. A. Vladimirov, Theor. Math. Phys. 43 (1980), 417. 

K.G. Chetyrkin, A.L. Kataev k F.V. Tkachov, Nucl. Phys. B174 (1980), 345. 

G. Curci k E. d'Emilio, Phys. Lett. B83 (1979), 199. 

M. Misiak k M. Miinz, Phys. Lett. B344 (1995), 308. 

K.G. Chetyrkin, M. Misiak k M. Miinz, Nucl. Phys. B518 (1998), 473. 

C. Ford, I. Jack k D.R.T. Jones, Nucl. Phys. B387 (1992), 373. 



J. A.M. Vermaseren, |math-ph/0010025 



P. Nogueira, J. Comput. Phys. 105 (1993), 279. 

S.A. Larin k J.A.M. Vermaseren, Phys. Lett. B303 (1993), 334. 

D.J. Gross k F.J. Wilczek, Phys. Rev. Lett. 30 (1973), 1343; H.D. Politzer, Phys. Rev. 
Lett. 30 (1973), 1346. 

D.R.T. Jones, Nucl. Phys. B75 (1974), 531; W.E. Caswell, Phys. Rev. Lett. 33 (1974), 244. 

O.V. Tarasov k A.A. Vladimirov, Sov. J. Nucl. Phys. 25 (1977), 585; E. Egorian k O.V. 
Tarasov, Theor. Math. Phys. 41 (1979), 863. 

L. Baulieu k J. Thierry-Mieg, Nucl. Phys. B197 (1982), 477. 

R. Alkofer k L. von Smekal, Phys. Rept. 353 (2001), 281. 

D.V. Nanopoulos k D.A. Ross, Nucl. Phys. B157 (1979), 273. 

R. Tarrach, Nucl. Phys. B183 (1981), 384; O. Nachtmann k W. Wetzel, Nucl. Phys. B187 
(1981), 333. 

K.-I. Kondo, T. Murakami, T. Shinohara k T. Imai, Phys. Rev. D65 (2002), 085034. 
K.G. Chetyrkin, S. Narison k V.I. Zakharov, Nucl. Phys. B550 (1999), 353. 



8 



